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Mixing properties of discrete-time quantum walks on two-dimensional grids with 
torus-like boundary conditions are analyzed, focusing on their connection to the 
complexity of the corresponding abstract search algorithm. In particular, an exact 
expression for the stationary distribution of the coherent walk over odd-sided lattices 
is obtained after solving the eigenproblem for the evolution operator for this particu- 
lar graph. The limiting distribution and mixing time of a quantum walk with a coin 
operator modified as in the abstract search algorithm are obtained numerically. On 
the basis of these results, the relation between the mixing time of the modified walk 
and the running time of the corresponding abstract search algorithm is discussed. 
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I. INTRODUCTION 

The concept of quantum walk is analogous to that of classical random walks, with the 
walker replaced by a quantum particle and coherent superpositions playing a key role 
There are discrete-time and a continuous-time versions of the quantum walk — the former 
defined by Aharonov et al. [2] and the latter, by Farhi and Gutmann jsj. Similarly to 
random walks, which have been used as a basis for classical algorithms that outperform their 
deterministic counterparts 0], quantum walks have also been used as a basis for quantum 
algorithms that outperform their classical correspondents j5|-[8[. 

Aharonov et al [9| have presented important results on the theory of quantum walks 
on graphs, several of them concerning mixing-time properties. In particular, they have 
analyzed odd-sided iV-cycles and showed that the quantum walk converges to a stationary 
distribution in time O(nlogn), almost quadratically faster than the classical walk. They 
have also obtained bounds on mixing-times for general graphs. Mixing times on hypercubes 
were properly addressed by Marquezino et al [10], who presented an analytical expression 
for the stationary distribution of a coherent discrete-time quantum walk on the hypercube, a 
Cayley graph with important algorithmic applications [5]. The stationary distribution is the 
first step in the calculation of the mixing time. An important graph is missing in the mixing- 
time picture: the two-dimensional grid. This grid, with torus-like boundary conditions, is 
the natural extension of the A^-cycle to two dimensions. The quantum search in this graph 



is one of the first examples of abstract search algorithms [11|, Il2| , a general framework for 
developing and analyzing quantum walk searches on graphs. 

In this paper, we consider the mixing time of a discrete-time quantum walk on torus- 
like 2D-grids, and analyze the relation between these properties and the complexity of the 
corresponding abstract search. First, the eigenvector problem of the evolution operator for 
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the quantum walk on the two-dimensional torus is solved. This is an essential result for 
the analysis of many mathematical properties of the quantum walk. Then, the limiting 
probability distribution of the coherent walk in odd-sided lattices is derived, for an initial 
condition localized at the origin in a uniform superposition of coin states. In this case, the 
limiting probability distribution is found to have a maximum at the origin. The mixing time 
has been numerically calculated and we show that the walk mixes in time 0{y/N log N), 
where N is the total number of vertices. 

A quantum walk in which the coin is modified according to the prescription of the abstract 
search algorithm has also been considered. Its limiting distribution and its mixing time 
have been numerically obtained. These results imply a relation between the mixing time 
of the modified walk and the complexity of the corresponding abstract search algorithm. 
In previous works, the complexity of the quantum search algorithm was considered to be 
related only to the hitting time, while the mixing time was considered to be related to data 
sampling 3| • We show how previous knowledge of the limiting distribution of the quantum 
walk may be used to estimate the running time of an abstract search algorithm. 

The paper is organized as follows. In Section [Tll we consider the coherent quantum 
walk on the two-dimensional finite grid with torus-like boundary conditions and solve the 
eigenvector problem of its evolution operator. In Section III Bl the limiting distribution of 
the coherent walk for the case of odd-sided lattices is derived. In Section IIII At the results 
of numerical simulations used to estimate the mixing time for the coherent evolution are 
presented. In Section IIII B\ the relation between the mixing time and the running time of 
the abstract search algorithm on the two-dimensional grid is discussed. In Section HVt the 
main results are summarized and our conclusions are presented. 



II. COHERENT WALK ON A TWO-DIMENSIONAL GRID 



A coined quantum walk in a viV x yiV grid with periodical boundary conditions has 
a Hilbert space Tic ^T-Lp, where He is the 4-dimensional coin subspace and Tip the A'"- 
dimensional position subspace. A basis for Tic is the set s)} for < rf, s < 1 and Tip is 
spanned by the set {\x,y)} with < x,y < a/A, where we assume integer \/N. A generic 
state of the quantum walk is 

1 v^-l 

1"^^)= E E ^d,s;.,ym,s)\x,y). (1) 
d,s=0 x,y=0 

The evolution operator for one step of the walk is 

U = S-{C(S)I), (2) 
where / is the identity in "Hp, S is the shift operator 

1 VAf-l 
d,s=0 x,y=0 

and C is a unitary coin operation in 'He- Notice that the binary sum © in 5* inverts the 



direction, as required by the abstract search algorithm We shall conisder the Grover 
coin, C = G, given by 

G = 2\u){u\-I, (4) 
where \u) = | ^^^=0 M' ^) uniform superposition in T-Lc- 
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A. Eigenproblem for U 

The analysis of the problem is simplified in Fourier space. Since the two-dimensional grid 
with periodic boundary conditions is a Cayley graph of we use the Fourier transform 
on this group, which has a basis spanned by the N kets 



^ VN-l 



\x,y), 



(5) 



x,y=0 



where uj = . The components of the evolution operator in the Fourier space are 



For each k^, ky, we define a reduced evolution operator in the coin subspace given by 



Gd,s;d',s' = J~'^y(-^''ok.+Sdiky)Q 



d,s®l;d',s' 



(7) 



which is a 4 X 4 matrix that can be diagonalized. The eigenvectors of U are tensor products 
of the eigenvectors of G and \kx, ky). 

Let us now describe the eigenspectrum of G. li k^ — and ky — 0, the eigenvalue 
1 is three-fold degenerate and the eigenvectors are — 1, 0, 0)^, ;^(0, 0, 1, — 1)^ and 

\u) = 1, 1, 1)^, where (. . .)^ is a column vector. The eigenvector with eigenvalue —1 is 
1, —1, —1)^. If /Ca; 7^ or ky 7^ 0, the eigenvalues are ±1 and e^'^^ with 6 defined by 



cos 9 — - 
2 



cos 



-|- cos 



27iL 



(8) 



In this case, the eigenvectors of G with eigenvalue +1 are 



4sin(^/2) 



1 — uj^y 



'■y (l — cu''^) 
uT'' - 1 



(9) 



and those with eigenvalue —1 are 



1 



I k^,ky/ 4cos(^/2) 
Finally, the eigenvectors with eigenvalues e*^ are 



-u;*^^ (l +uj'^y) 
- (1 +0;*^^) 
u''y (1 + w''^) 
1 + 0;^=- 



(10) 



2^2 sin 61 



e~*^ — cj'^^ 

e"*^ — uj^y 
e~*^ — u"^'" 



(11) 
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and the eigenvectors with eigenvalue e are obtained by the replacing 6 — 5- —6 in Eq. (Illl) . 
Note that the eigenvectors are normalized, \i^k^ky)i fo^^^^ ^in orthonormal basis for the reduced 

space and have a real constant component on the uniform state, {y'^^^J\u) = 1/a/2. 
We take the state 

|*(0)) = |n)|a; = 0,y = 0) (12) 

as the initial condition, i.e., the walker starts localized at the point (0,0) and uniformly 
distributed in the coin subspace. In the eigenbasis, this initial condition is given by 

1 1 

1^(0)) = ^ \u)\K = 0,ky = 0) + -^ Yl + l^-^J/)- (13) 

Applying [/* on |^E'(0)) we obtain the state of the quantum walk after t steps, 
1 1 

(A;^,fcy)^(0,0) 

Having solved the eigenvector problem for the evolution operator, we proceed to the calcu- 
lation of the limiting distribution for the quantum walk. 



B. Limiting Distribution 

Let P{x, y, t) be the probability to find the walker at a vertex {x, y) of the grid at time t. 
As mentioned in the introduction, this probability depends on the initial condition and, as 
is typical of unitary evolutions, it does not converge to a stationary distribution. However, 
the time-averaged distribution P{x, y,T) = ^ Ylt=o P{^^ ^) always converges as T goes to 
infinity [9]. The limiting or stationary distribution is then defined in terms of the average 
distribution, 

7i"(a;, y) = lim P{x, y, T). (15) 



Using Theorem 3.4 from Ref. [9|, the coefficients of the initial condition expressed in the 



.... \.'+^' 1"+^ 

expression for the limiting distribution 



eigenbasis of U, Eq. (IT^ . and the fact that (z/^ fc' I'^aL fc^) ~ ('^k' y I'^k^ ky)i obtain a simple 



iV2 iV2 

(fc:r,fcy)7^(0,0) 

(fc;,fc;)^(o,o) 
e{k^,kv)=e{k'^,k'y) 



When 6{k'^, k' ) = 9{kx, ky) we have 



I 6 I e \ _ 1 - 2 cos^ [9{kx, ky)] + cos [O^k^ - k'^, ky - k'y]^ 
{^k'^,k'J^k.,ky) - 2sm' [e{kx,ky)] • 
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Analyzing all cases such that 0{k'^, ky) = 6{kx, ky), the expression for the limiting distri- 
bution, for odd V^, has the explicit form 

,1 2 1 

X I [(x - y) cos k, + co'^-^'^+y^j (l + cos Ai,) + 2 (w^^.- + ^^k^y^ | + 



2A^2 



VN-l 

E 



COS {kx — ky) — cos2^^ 



+ 



cos^ kx — cos 26 



'2kxX 



cos ky — COS 26* 



(18) 



+ [cos — ky, kx + fc^)) — cos2^^] u 
1 



+ - 



sin kj. + sin A;^ 



1 
2 

+ 



cos A;^; — cos ky 



(kx-ky)x+(kx+ky)y 
{kx~\~ky)x-\-{hy kx^y 



COS {kx + ky) — cos 29 



UJ 



{kx~\'ky^x~\-{^kx~\~ky^y 



where kx = 2'Kkx/vN and ky = 2'Kky/ y N . 

For some values of (x, y), it is possible to simplify Eq. (IT5]) and achieve simple results. An 
interesting point to be considered is the initial site, which was assumed to be (xq, yo) = (0, 0), 
without loss of generality. It is straightforward to show that 



7r(0,0) 



AN -8VN + 5 
iV2 



(19) 



and, for ^ 1, we obtain vr(0,0) ~ A/N, which is also the maximum of the limiting 
distribution. This distribution is shown in Fig. [T]for a grid of dimensions 41 x 41, with the 
initial site shifted to {xo,yo) = (20,20) for better visualization. 

It is well known that the behavior of quantum walks on even lattices may be different than 
that observed for odd lattices jo], [l3]- The analysis of the walk over odd lattices is sufficient 
for the main objective of this paper, namely to study the relation between mixing time and 
the complexity of quantum-walk based search algorithms. Our numerical simulations show 
that the liming distribution for even lattices present two peaks, instead of the single peak on 
the initial site as observed for the odd lattice. Now that we have the limiting distribution for 
the quantum walk on the two-dimensional grid, we consider its mixing time and its relation 
to the abstract search problem. 

III. MIXING TIMES AND ABSTRACT SEARCH 



In this subsection we consider the mixing time for a coherent evolution. The rate at 
which the average probability distribution of a quantum walk approaches its asymptotic 
distribution is captured by the following definition Q. 
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Figure 1: Left panel: Limiting distribution for quantum walk in two-dimensional grid with ^/ N = 
41, obtained from Eq. ()18p with the initial condition (|12p . Right panel: contour plot for the same 
distribution. 



Definition III.l The average mixing time of a quantum Markov chain to a reference 
distribution vr is 

M, = mm{T | Vt > T, \\Pt - n\\ < e}, 

where \\A — B\\ = \A{x) — B{x)\ is the total variation distance between the two distribu- 
tions. 

An alternative definition captures the first instant in which the walk is e-close to the reference 
distribution tt, 

Definition III. 2 The instantaneous mixing time 1^ of a quantum Markov chain is 

= min{t I \\Pt — 7i\\ < e}. 

Both mixing times depend on the initial condition of the quantum walk. 

For a quantum walk in a generic graph with arbitrary initial condition, an upper bound for 
the total variation distance to the asymptotic distribution 7i{x,y) was derived by Aharonov 
et al [9|, 

TT 

\\P{x,y,T) - Tx{x,y)\\ < 



TA 



log ( — + 1 



(20) 



where N is the number of vertices, d is the degree of each vertex and A is the minimum 
separation between distinct eigenvalues of U. 



A. Mixing time in the two-dimensional grid 

Let us now consider the particular case of a two-dimensional cartesian grid, for which 
d = A. The value of A is the minimum value of \^e^^'^^^''^v) — Q^^C^^^ky)^ -^^ 

the range [O, — l] . For large and small values of k^, ky, k'^ and ky, we obtain 

V2n 



A 



^2 I 1.2 



1.12 \ y2 



(21) 
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The minimum value of A is obtained taking values of k^, ky, k'^ and k'y such that 9{kx-, i^y 
OiK, K) and 



1.2 _|_ 1.2 



U2 \ U2 



We can make a good guess by exploring the structure of the above equation. With those 
guesses we can try to find an upper bound for the average mixing time, M^. It is easy 
to obtain the value 0(1/V^) for A. Therefore, a good guess for is O ^y^El^sK^^ gy 

performing a numerical analysis, we have succeeded in obtaining better bounds regarding the 
dependence on N. Plotting against ^/TT\ogN for several values of e we obtain straight 
lines as shown in the right panel of Fig. [2J These results show that is proportional to 
where c is approximately 1. Therefore, the numerical data strongly suggests that 



M, = e 



^N\ogN 



(22) 



In the left panel of Fig. |2l we have the total variation distance to both the uniform and the 
stationary distribution. For long times, the variation distance to the stationary distribution 
decays approximately as ~ 1/t while the corresponding distance to the uniform distribution 
remains essentially constant. 
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Figure 2: Left panel: total variation distance to both the uniform and the stationary distributions 
of the quantum walk on the two-dimensional grid with flip-flop shift as a function of the time step. 
Right panel: mixing time to the stationary distribution as a function of the input size. 

The classical random walk on a two-dimensional grid increases with the size of the lattice 
as Q{N) — see, for instance, chapters 3 and 5 of Aldous and Fill [15|]. Comparing this 
result with Eq. f l22|) . we observe that the quantum walk mixes almost quadratically faster 
than its classical counterpart on the same lattice. The faster mixing rates observed in the 
quantum case is one of the main advantages of using quantum algorithms over their classical 
equivalents. 

Until now, we have analyzed the standard quantum walk on the torus, without direct 
algorithmic applications. In the next section, we investigate the behavior of a quantum walk 
with a modified evolution operator which makes it useful to to mark a searched vertex of 
the grid. 
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B. Connection with abstract search algorithms 

In this section we discuss the relation between the mixing time and the running time 
of the abstract search algorithm on the two-dimensional grid. The running time of search 



algorithms is usually associated with the notion of hitting time [16[, which is defined as the 
first time a given vertex is reached. On the other hand, the mixing time is usually associated 
with data sampling [l3[. We argue that those measures are related and both may be used 
to estimate the running time of abstract search algorithms. 



The abstract search algorithm is a search framework introduced by AKR based on 
a modified quantum walk. The standard quantum walk is driven by the evolution operator 
given by Eq. The modified operator is U' = S ■ C, where C is given by 

C" = -J® \xo,yo){xo,yo\ + G (g) (/ - \xo,yo){xo,yo\). (23) 

The modified operator C applies the coin — / if the vertex is the target |xo, yo), otherwise it 
applies Grover's coin operation G. The effect of C is to mark the searched vertex with a rel- 
ative phase. With this new evolution operator, AKR have shown that a quantum walker de- 
parting from the uniform distribution will be at the marked vertex after 0{y/N log N) steps 
with probability 0(1/ log A^). The time complexit y o f AKR's algorithm is 0{\/N\ogN), 
after using the method of amplitude amplification |l7l |. Tulsi [l2| has improved this result 
and have shown a method to reach the marked vertex with probability 0(1). The time 
complexity of Tulsi's algorithm is 0{y/N log A^). 

The effect of the modified evolution operator is to increase the probability of finding the 
walker in the marked vertex at some specific steps. The best point to stop is the one at which 
the probability is maximum. AKR described a method to find the running time without 
calculating directly the point of maximum and without using the notion of hitting time. 
On the other hand, Szegedy [l6| has described a general method to obtain the hitting time. 
Szegedy's method, however, requires a change on the evolution operator, which does not have 
the same time complexity of AKR's method for searching one vertex on the two-dimensional 
grid. 
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Figure 3: Quantum walk in two-dimensional grid with ^/N = 41 and a modified coin used to 
search for a marked node. Left panel: Probability distribution after t = 80 steps, corresponding 
to the instant of maximum probabiUty at the marked node. Right panel: stationary distribution 
approximated with T = 10^ simulation steps. 

The left panel of Fig. [3] displays the probability distribution for the case \fN = 41 after 
80 steps, which corresponds to the first maximum of the probability at the marked vertex. 
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The running time is the instantaneous mixing time using as reference distribution the one 
with maximum probabihty at the marked vertex and taking a small value of e. This fact 
cannot be used to estimate the running time, because there are no clues about how to find 
that specific probability distribution beforehand. On the other hand, there is a method to 
find a priori the limiting distribution. Hence, an interesting question is whether the limiting 
distribution can be used to estimate the running time. The right panel of Fig. [3] displays 
the stationary distribution. Note the remarkable similarity between the distribution with 
maximum probability and the stationary distribution. 

The first thing we need to check is how the mixing time of the modified walk scales with 
A^. The right panel of Fig. H] strongly suggests that replacing C by C defined in Eq. fl25]) 
does not alter the scaling of the average mixing time, since for ^ 1 it scales as in eq. fl2^ . 
There is a small oscillation in the mixing time function for large N, related to the way the 
modified walk approaches the limiting distribution. In the the left panel of Fig. IH we have 
the total variation distance to both the uniform and the stationary distribution. For long 
times, the variation distance approaches the stationary distribution with a strong oscillation 
and several local minima, corresponding to the steps when the instantaneous probability 
distribution has a maximum at the searched node. 




Figure 4: Left panel: total variation distance to both the uniform and the stationary distributions 
of the quantum walk in two-dimensional grid with modified coin used to search for a marked node. 
Right panel: mixing time to the stationary distribution as a function of the input size. 



This behavior is remarkably different from the one observed in the last section for the 
Grover walk on the two-dimensional grid (Fig. |2]). In that case the mixing time increases 
without oscillations and the total variation distance decays almost as a power law. 



IV. DISCUSSION 



In this paper, the mixing properties of a discrete-time quantum walk on the finite two- 
dimensional grid with torus-like boundary conditions has been considered in detail. This 
particular topology is quite interesting due to its algorithmic applications. The relation 
between mixing time and the time complexity of the corresponding abstract search algorithm 
has also been investigated. 

The eigenvector problem of the evolution operator for the coherent quantum walk has 
been solved and the corresponding stationary distribution has been found analytically for the 
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particular case of odd lattices and a localized initial condition with a uniform superposition of 
coin states. The stationary distribution is not uniform and has a maximum at the origin — for 
which we have also provided a simplified expression. According to our numerical simulations, 

the mixing time on the two-dimensional grid with vertices increases as O ' ^ 



i.e., almost quadratically faster than the mixing time of the classical random walk, which is 
Q[N). The value of c is approximately 1. The faster mixing rates on the quantum case is 
an advantage of using quantum-walk based algorithms over the classical equivalents. 

We have also considered a quantum walk in which the coin is modified according to the 
prescription of the abstract search algorithm. We have numerically calculated its limiting 
distribution and its mixing time. Our numerical simulations show that the mixing time 

for this particular walk also increases as O(^^^l^), where c is approximately 1. This 

mixing time of the modified quantum walk corresponds to the time complexity of marking 
a searched vertex on the torus using the abstract search algorithm jl2|. This result estab- 
lishes a relation between the mixing time of the modified walk and the running time of the 
corresponding abstract search algorithm — which previously has been associated only to the 
hitting time. This relation may be useful for the complexity analysis of quantum-walk based 
algorithms. 

In order to obtain the exact running time of the quantum algorithm, one should also 
consider the constants that have been absorbed by the asymptotic notation used to represent 
the mixing time. As a future work, it would be interesting to investigate the mixing time of 
search algorithms on other graphs and compare the corresponding constants. 



Acknowledgments 

The authors acknowledge financial support from CNPq (Brazil) and PEDECIBA 
(Uruguay). 



[1] J. Kempe, Contemp. Phys. 44, 307 (2003), 'arXiv:quant-ph/0303081vl. 
[2] Y. Aharonov, L. Davidovich, and N. Zagury, Phys. Rev. A 48, 1687 (1993). 
[3] E. Farhi and S. Gutmann, Phys. Rev. A 58, 915 (1998), arXiv:quant-ph/9706062V2. 
[4] R. Motwani and P. Raghavan, Randomized algorithms (Cambridge University Press, UK, 
1995). 

[5] N. Shenvi, J. Kempe, and K. B. Whaley, Phys. Rev. A 67, 052307 (2003), 

arXiv:quant-ph/0210064 , 
[6] A. Ambainis, Quantum walk algorithm for element distinctness, in Proceedings 45th Annual 

IEEE Symp. on Foundations of Computer Science (FOCS), 2004, ar Xiv : quant-ph /0311001. 



[7] A. Childs and J. Goldstone, Phys. Rev. A 70, 022314 (2004), arXiv:quant-ph/0306054 



[8] E. Farhi, J. Goldstone, and S. Gutmann, A quantum algorithm for the hamiltonian NAND 

tree, arXiv:quant-ph/0702144v2. 
[9] D. Aharonov, A. Ambainis, J. Kempe, and U. Vazirani, Quantum walks on graphs, in 

Proceedings of 33th ACM Symposium on Theory of Computation (STOC'Ol), pp. 50-59, New 



York, NY, 2001, ACM, |aBCw:quant-ph/0012090| y2. 



11 



[10] F. Marquezino, R. Portugal, G. Abal, and R. Donangelo, Phys. Rev. A 77, 042312 (2008), 
karXiv:0712.0625. 

[11] A. Ambainis, J. Kempe, and A. Rivosh, Coins make quantum walks faster, in Proceedings of 
the sixteenth annual ACM-SIAM symposium on discrete algorithms, pp. 1099-1108, Society 
for Industrial and Applied Mathematics Philadelphia, PA, USA, 2005. 

[12] A. Tulsi, Phys. Rev. A 78, 012310 (2008), 0801.0497. 

[13] F. Magniez, A. Nayak, P. Richter, and M. Santha, On the hitting times of quantum versus 
random walks, in Proceedings of the twentieth Annual ACM-SIAM Symposium on Discrete 
Algorithms, pp. 86-95, New York, NY, 2009, Society for Industrial and Applied Mathematics. 

[14] V. Kendon and B. Tregenna, Phys. Rev. A 67, 042315 (2003), |arXiv :quant-ph/0209005[ 



[15] D. Aldous and J. Fill, Reversible Markov Chains and Random Walks on Graphs, 

jhttp:/ / stat-www.berkeley.edu/users / aldous /RWG7book.html , 1995. 
[16] M. Szegedy, Quantum speed-up of markov chain based algorithms, in Proceedings of the 45th 

Annual IEEE Symposium on Foundations of Computer Science, pp. 32-41, 2004. 
[17] P. Kaye, R. Laflamme, and M. Mosca, An Introduction to Quantum Computing (Oxford 

University Press, Inc., New York, NY, USA, 2007). 



